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^ ' Abstract 

^ ' We study the mean-field dynamics of a fermionic condensate interacting with a single bosonic 

. mode (a generalized Dicke model). This problem is integrable and can be mapped onto a corre- 

sponding BCS problem. We derive the general solution and a full set of integrals of motion for 
C ' the time evolution of coupled Fermi-Bose condensates. The present paper complements our earlier 

^ ■ study of the dynamics of the BCS model (cond-mat/0407501 and c ond-mat / 0505493 1 ) . Here we 

provide a self-contained introduction to the variable separation method, which enables a complete 

> 

CN ' analytical description of the evolution of the generalized Dicke, BCS, and other similar models. 
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I. INTRODUCTION 



A number of closely related non- stationary problems have come up recently in different 
contexts. In these problems the goal is to describe the dynamics of a many-body system 
following a sudden perturbation that drove the system out of an equilibrium. The system 



in question can be a BCS superconductor 



Bose condensates 







ll|, coupled Fermi- 



13 L o r a sin gle electronic spin interacting with many nuclear spins 
(the central spin model) ij, IigIIitI lisl liol Q • A common feature of all these problems 
is that they can be formulated in terms of spin Hamiltonians, which belong to a class of 



integrable systems known as Gaudin magnets 



21 



22 



3- 



It turns out that this fact enables 



one to solve for the evolution of all these systems using the same approach. 

In the present paper we focus on a particular example of such a problem. Namely, 
we study the dynamics of a fermionic condensate interacting with a single bosonic mode. 



Recent studies of this problem 



12 



llS^ were motivated by exDerimei itsl24l. 

HQ 



26 



27| on 



cold fermion pairing. The idea was that near the Feshbach resonance j2Sl, |22| these systems 
can be modeled by a fermionic condensate of atoms strongly coupled to the lowest energy 
bosonic mode. 

We assume that the system has been prepared in a nonequilibrium state at t < and 
study the subsequent evolution for t > 0. Our main result is a complete solution for the 
dynamics of the system. We also derive a full set of its integrals of motion. It turns out that 
the mean-field time evolution of the coupled Fermi-Bose condensates can be mapped onto 
the corresponding BCS problem, which was solved previously in Refs. For this reason, 

the emphasis in the present paper is on this mapping and detailed or alternative derivations 

9, 



that have been largely omitted in Refs. 



IG 



The fermion-boson condensate is described by the following Hamiltonian: 



(1.1) 



where are the single-particle energy levels and the operators c]^. (cjo-) create (annihilate) a 
fe.n,.on of one of tl>e two species . =T o. i in an o.bital eigenjtate of ene.gy Egen.tate. 
\j t) and \j I) are related by the time-reversal symmetry j30|]. For example, if the single- 
particle potential is translationally invariant, |j ]) = |p |) and |j i) = | — p i)- Operators 
(6) create (annihilate) quanta of the bosonic field. 
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We study the dynamics of the fermion-oscillator model in the mean-field approxima- 
tion. This amounts to treating the bosonic field classically, i.e. replacing operators and b 
with c-numbers in the Heisenberg equations of motion for Hamiltonian (jl.lj) . This procedure 
is expected to be exact as long as the bosonic mode is macroscopically populated. It turns 
out that in this approximation the dynamics coincides with that of a classical Hamiltonian 
system, which can be mapped onto the corresponding BCS problem (see below). The clas- 
sical dynamical variables are the time-dependent quantum-mechanical expectation values 
(cj|Cj|), (b), and (Z^o- c]^Cjo-). We will see below that the mean-field approximation is also 
equivalent to the usual way of obtaining the classical limit - by replacing operators with 
classical variables and their commutators with Poisson brackets. 

Here we assume that the number of energy levels, n, interacting with the bosonic field 
in the system (jl.ll) is arbitrary large, but finite. In this case the typical evolution is 
quasi-pehodic with n + 1 incommensurate frequencies. The system exhibits irregular multi- 
frequency oscillations ergodically exploring the part of the phase-space allowed by the con- 
servation laws, returning arbitrarily close to its initial state at irregular time intervals. For 
finite n, the solution, albeit explicit, is rather complicated. However, it considerably sim- 
plifies in the thermodynamic limit n ^ oo. In this limit, the return time diverges for 
most physical initial conditions, while the solution decays to a simple limiting dynamics 
at large timesjs^. For exarnple, following an initial decay, \(b{t))\ asymptotes either to a 



constant value (cf. Refs. y, |3|, lH, [11 ) or to an oscillatory behavior characterized by only a 
few frequencies (cf. Fig. 2 in Refs^,y). 

The fermion-oscillator model (jl.lj) can be viewed as a generalization of the Dicke (Tavis- 
Cummings) model of Quantum Optics. The latter can be obtained from (jl.lj) in the zero 
fermionic bandwidth limit, i.e. when all single-particle levels ej are degenerate, ej = fi. To 
see this, it is useful to reformulate the model (jl.lj) as a spin-oscillator model using Anderson's 

. Pseudospins are defined as 

2k^ = n, - 1, k- = c,^c,^, = c]^c]^, (1.2) 

where Uj = Y.a^^jaCja- Operators have all properties of spin-1/2 on the subspace of 
unoccupied and doubly occupied (unblocked) levels e^. Singly occupied (blocked) levels 
do not interact with the bosonic field and are decoupled from the dynamics. In terms of 
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pseudospins the Hamiltonian p.ip for n unblocked levels takes the form 

n— 1 n— 1 

H=Y. 2e,^| + ojh^h + gY, [pK' + bK+) . (1.3) 

j=0 j=0 

In the zero bandwidth limit the Hamiltonian (|1.3|) describes an interaction of a single 
collective spin T = J2j K-, with a harmonic oscillator 

Hoicke = 2fit + ujPb + 9{b^f_ + bf+). (1.4) 

This model was suggested by R. H. Dicke in 1953 js^. Its spectrum was obtained exactly 
by M. Tavis and F. W. Cummings in 1967 33j. In the mean- field approximation it becomes 
a classical model for the dynamical variables (T(t)) and {b{t)) with only two degrees of 
freedom (see also below). The mean-field solution for the time evolution was outlined by 

n n 

Dicke j32j and derived in detail by R. Bonifacio and G. Preparata in 1969 j3j]. In fact, the 
resulting classical problem is that of a spherical pendulum. The solution for (T(t)) and (bit)) 
is in terms of elliptic functions and (apart from the azimuthal motion of the pendulum) is 
simply periodic |34| . 

Remarkably, a more general many-body model ()l.lll.3p also turns out to be integrable. 
This was established by M. Gaudin in 1972 2^ and later by B. Jurcojs^ who used a dif- 
ferent approach. First explicit nonlinear solutions for the mean-field dynamics of the model 
()l.lll.3j) were constructed in Refs. Interestingly, the initial conditions for which these 

solutions occur are such that the dynamics of the model ()l.lll.3|) reduces to that of the Dicke 
model ()1.4j) . i.e. it can be described in terms of a single collective spin coupled to the bosonic 
field. Similar solutions for the mean-field dynamics of the BCS model were discovered by 
R. A. Barankov et. al. in Ref. and in an unpublished work by V. S. Shumeiko P]. 

Here we solve for the dynamics of the model (jl.lll.3p in the mean-field approximation for 
arbitrary initial conditions. We also derive a full set of conservation laws for the mean-field 
dynamics. The typical evolution of the system is quasi-periodic with n + 1 incommensurate 
frequencies. However, for certain special choices of initial conditions the dynamics is char- 
acterized bym-|-l<n + l incommensurate frequencies and can be described in terms of m 
collective spins coupled to the bosonic field. 

Our approach to the mean-field dynamics of the fermion-oscillator model (jl.lj) employs 
the method of separation of variables. This method was suggested by I. V. Komarovj37| and 



later developed by E. Sklyanin and V. Kuznetsov 



22, It allows us to derive a full set of 



integrals of motion for both the quantum model ()l.lll.3p and its classical counterpart. We 
also use it to derive and analyze equations of motion in terms of separation variables. 

Upon the replacement of the quantum bosonic field b with its time-dependent expectation 
value, the Hamiltonian becomes similar to the mean-field BCS Hamiltonian. In this 
analogy, g{b{t)) plays the role of the BCS gap A(t). The difference is that in the case of the 
BCS model A(t) is not an independent dynamical variable, but is related to {cji{t)cj^(t)) = 
{Kj{t)) by the self-consistency condition A(t) = gJ2j{K^(t)). Nevertheless, it turns out 
to be possible to map the mean-field dynamics of the fermion-oscillator model p.lj) with 
n energy levels onto the corresponding BCS problem with n + 1 levels. This mapping 
is facilitated by the variable separation and enables us to obtain the general solution for 
the time dependence of the expectation values {Kj{t)) and {b{t)) in terms of hyperelliptic 



39| using the known solution for the mean-field dynamics of the BCS model 



functions 38. 

Q. 

The rest of the paper is organized as follows. In section |nj we derive the classical 
Hamiltonian that governs the mean-field dynamics of the quantum model ()l.ipi.3j) . We 
show that the dynamics of the resulting classical model is integrable and derive a full set 
of its conservation laws in section IIIII In section IIVI we perform a transformation to a new 
set of variables, which facilitates the solution of the equations of motion. The mapping to 
the corresponding BCS problem and the general solution for the mean-field dynamics of the 
model ()l.lll.3|) are derived in section IVIl Finally, section IVIII is devoted to an interesting 
class of particular solutions that include mean-field equilibrium states and special solutions 
of Refs. Ill Q, HQ- 



II. CLASSICAL MODEL 

Here we derive the classical Hamiltonian model that describes the mean-field dynamics 
of the quantum model ()l.lll.3|) . 

We start with the Heisenberg equations of motion for the spin-oscillator model ()1.3|) . 



Kj = i H, Kj = Bj X Kj Bj = {2gb^, 2gby, 2e 



n-l 



(2.1) 



b = t[H,b] = -lujb - igT_ T = 5^ Kg, 

q=0 
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where the operators bx and by are defined by 

b = bx — iby b^ = bx + iby. 

In the regime when the bosonic mode is macroscopically populated, we can replace operators 
b{t) and 6^(t) in Eqs. ^2,.\\ with c-numbers: b{t) = (b{t)) and b{t) = (b^{t)), where (. . .) stands 
for the time-dependent quantum- mechanical expectation value. After this replacement, 
Eqs. ()2.H1 become linear in operators. Taking their quantum-mechanical expectation value, 
we obtain for Sj{t) = (Kjit)) 

Sj = Bj X Sj Bj = {2gbx, 2gby, 2tj), 

(2.2) 

'i- 



b = —iujb — igJ^ J ~ 51 

j=0 



An important observation is that Eqs. ()2.2|) are Hamiltonian equations for the following 
classical model: 

n—l n—l 

H=Y^ 2e,s] + ubb + gJ2 (bsj + fo+) , (2.3) 

j=0 j=0 

where sf = ±^sj. This Hamiltonian governs n classical spins interacting with a classical 
harmonic oscillator j^. 

Classical spins are similar to the angular momentum and have the same Poisson brackets 

(2.4) 

{b,b} = i {bx,by} = -, 

where a, b, and c stand for the spatial indexes x, y, and z. All other Poisson brackets 
between components of and b vanish. Equations (|2.4j) determine the fundamental Poisson 
brackets in our problem in a sense that Poisson brackets between any other pair of dynamical 
variables (functions of b,b and components of s^) can be obtained from Eqs. ()2.4|) using the 



standard properties of Poisson brackets 4l(]. For example, Eqs. fj2.2|) can be obtained from 
the Hamiltonian equations of motion Sj = {H, Sj} and b = {H, b} using Eqs. ()2.4|1 and ()2.3|1 . 

Equations of motion (|2.2p conserve the length of each spin, s| = const. This fact can also 
be viewed as a property of the brackets (j2.4j) - since Poisson-commutes with all other 
dynamical variables, its bracket with the Hamiltonian ()2.3|) also vanishes. In the mean-field 
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approximation eigenstates of the Hamiltonian p.3|) are product states. In a pure state of a 
spin-1/2, there is always an axis, n, such that the projection of the spin onto it is 1/2, i.e. 
Kj ■ n = 1/2. In this case, s| = (K^)^ = 1/4. Therefore, if the system was in a product 

state at t = 

= i^'^r + stsj = \. (2.5) 

If a number of orbitals ej in ()2.3|) are degenerate, the magnitude of their total classical 
spin J2ej=e is conserved by the equations of motion ()2.2|) . In this case, one can replace the 
corresponding spins Sj with a single classical spin of a larger magnitude, 

ejSj = eSe, where = ^ Sj, (2.6) 

in the Hamiltonian ()2.3|) and sum over nondegenerate orbitals only. Below we will assume 
whenever necessary that such a replacement has been made, i.e. that orbitals ej are nonde- 
generate. 

Finally, comparing Hamiltonians ()1.3p and ()2.3p and brackets fl2.4j) to the corresponding 
quantum commutators, we note that the mean-field approximation is equivalent to the 
standard procedure of going from quantum to classical mechanics by replacing commutators 
with Poisson brackets, i[v4,i?] {A,B} (here h= 1). 



III. INTEGRABILITY 

Here we demonstrate that the classical model ()2.3|) . as well as its quantum counterparts 
fll.3|l and p.lj) . are integrable and introduce a useful tool to analyze and solve for their 
dynamics. 

The Hamiltonian ()2.3|) depends on 2{n + 1) dynamical variables: 2 angles for each of n 
spins plus the coordinate and the momentum of the oscillator j^^. Therefore, its phase space 
is 2{n + l)-dimensional and the number of the degrees of freedom (the number of generalized 
coordinates) is n + 1. To show that the classical model ()2.3|) is inte grab le, we have to show 
that it has n + 1 independent integrals of motion (see e.g. Refs. 0, 111). 
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Consider the following vector-functions of an auxiliary parameter u: 

Si 



Lj{u) 



u — e,- 



j = 0, . . . , n - 1 



/ 2gh, \ 
2gby 
2u — uj 



It follows from the above definitions and Eqs. ()2.4j) that components of Lfe(u) have the 
following Poisson brackets: 



^abc 
V — W 



{Ll{v) - Ll{w)) k = 0,...,n 



(3.1) 



Because components of different Lj('u) Poisson-commute with each other, their sum 



Uu) 



2gby 
2u — UJ 



n-l 



+ E 



j=o « - 



(3.2) 



also satisfies Eqs. (j3.1|) . One can check using only Eq. (j3.ip that the square of this vector 
commutes with itself at any two values, v and w, of the auxiliary parameter, i.e. 



(3.3) 



The function L^('u) acts as a generating function for the model ()2.3|) and its integrals of 
motion. Indeed, evaluating L^(m) from Eq. ()3.2|) . we obtain 



L (w) = -i + —^ + 1^ 

9 ^9 j=o 



2Hi 



where 



Hj = {2ej - u)s] + gibs- + bs+) + g^Y. 



Sj • Sfc 



Hr, = UJ 



(bb + J, 



n-l 

i=o 



(3.4) 

(3.5) 
(3.6) 



Hamiltonians Hj defined by Eqs. ()3.5|3.(ij) Poisson-commute with each other 



{iJj, Hk] = for all j and k. 
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This follows directly from the fact that Eq. ()3.3p holds for any v and w. On the other 
hand, it is straightforward to verify that the spin-oscillator Hamiltonian ()2.3p is a linear 
combination of Hj 

n 

H = Y.Hr (3.7) 

Since Hj Poisson-commute with each other, the Hamiltonian ()2.3|) commutes with all n + 1 
Hamiltonians Hj. Therefore, Hj given by Eqs. ()3.5l3.6p are integrals of motion for the 
classical spin-oscillator model, i.e. this model is integrable. 

The above construction based on Eqs. (|3.ip and (|3.2|) constitutes the so-called Lax vector 
representation for the spin-oscillator model (j2.3j) . Its main advantages are that it provides 
a powerful tool to study the dynamics and can also be extended to quantum models. In 
particular, note that integrals ()3.5|) and ()3.6p can be trivially generalized to the quantum 
case 

H^ = (2e, - + g^Pkj + + g^Y. 



j^k 



n-l 



Hr, = u;(b^b + t) T = Y.K,. 

q=0 

Operators Hj pairwise commute with each other and the spin-oscillator Hamiltonian ()1.3p . 
which is their linear combination as in Eq. (|3.7|) . 

The same construction with classical dynamics in the mean-field approximation and Lax 
vector for resulting classical models applies to a number of other models. In particular, a 
closely related model is the BCS model 

Hbcs = ^j^l^cja - qY. c]^c]^CqiCq^. (3.8) 

i.o- j,q 

Notations here are the same as in Eq. In terms of Anderson's spins (see above 

Eq. ()1.3|) ). the Hamiltonian reads 

n 

Hbcs = E ^^i^j ' 9 E (3-9) 

j=0 j,q 

The usual BCS mean-field is equivalent to the procedure (Eqs. ()2.1|) and ()2.2|) ) we used to 



derive t 
Refs 



" I lie 



le classical Hamiltonian ()2.3|) . only now the role of b is played by J2j Kj (see e.g. 
The mean-field dynamics is described by the following classical Hamiltonian: 

n 

Hbcs = E 2e,5 ■ - qY (3-10) 
9 



with the Lax vector 

LBcsiu) = --+^-^, (3.11) 
g jU-ej 

where z is a unit vector along the ^-axis. The Lax vector Lbcs(w) has the same properties 
fl3.1|) and ()3.3|) as before and generates integrals of motion for the BCS model. 



IV. SEPARATION OF VARIABLES 

In this section we perform a transformation to a new set of variables, which facilitates 
the solution of the equations of motion ()2.2|1 . The choice of new variables naturally follows 
from the Lax vector construction of the previous section. The new variables are canonical 
and also separate the Hamilton- Jacobi equations for the classical Hamiltonian ()2.3|1 in the 



usual sense 



3, 



We define n variables Uj as zeros of 

2b ""^ 

L_{u)=L^{u)~iLy{u) = - + Y.^^. (4.1) 

Note that the nominator of this expression is a polynomial of degree n and therefore there 
are n roots Uj with j = 0, . . . , n — 1. The coefficients of this polynomial are functions of 
the dynamical variables sj and b. Accordingly, its roots are also functions of sj and b and 
therefore define a new set of dynamical variables. 

Variables Uj play a role of canonical coordinates for the classical Hamiltonian ()2.3|) . The 
corresponding momenta are defined as Vj = Lz{uj). Thus, we have 

L_(n,)=0, v.j=L,{uj), j = 0,...,n-l. (4.2) 

Because our system has 2(?2 + 1) degrees of freedom, we need two additional variables Un 
and Vni which can be introduced as 



Un = b, Vn = ^ = b+ ^, (4.3) 
Lob b 

where Hn is given by Eq. ()3.6|) . 

Separation variables {uj,Vj) are canonical, i.e. 

{uj,Uk}=0 {vj,Vk} = {vj,Uk} = -iSjk j,k = 0,...,n. (4.4) 
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The first relation in Eq. ()4.4|) follows from the fact that by Eqs. ()4.2j) and ()4.3p variables 
Uj depend only on mutually Poisson-commuting variables and b. The second relation 
follows from {Lz{v), Lz{w)} = {Lz{v),b} = 0. To derive Poisson brackets between Vj and 
Uk for j, k = 0, . . . ,n — 1; we use the following equation obtained from Eq. (j3.ip 

{L^iv), L_H} = (L^v) - L_{w)) . (4.5) 

Evaluating this expression at v = Uj and w = Uk, we obtain 

{L^{uj),L^{w)}uj=Uk = — z — L_{w) = iL'_{uk)5jk- (4.6) 
where L'{u) = dL/du. On the other hand, 

{L^{uj),L^{w)}.^=u^ = {Lz{uj),L^{uk)} - L'_{uk){Lz{uj),Uk} = -L'_{uk){L^{uj),Uk}, 

(4.7) 

where we used L_^{uk) = 0. Comparing Eqs. ()4.6|) and ()4.7|) . we obtain the last relation in 
Eq. dOl) for j. A; = 0, . . . , n - 1. 

The original djTiamical variables = (Kj(t)) and b = {b(t)) that we are interested in, 
can be explicitly expressed in terms of uj (the inverse map). Indeed, note that L^{u) has 
its zeros at m = Uj and poles at m = tj. Therefore, using Eqs. ()4.2|1 and (jHH), we can write 
it as 

9 ,t; " - 9 nj(« - «i) 

where we have also used L_{u) = 2b/g + 0{l/u) for large u. Equating residues at m = 
and M = oo, we obtain 

= - -^-7 (4.9) 



2& nfc(e,- -»fc) 



2b 

Similarly, using 

L^{uj)=Vj L^{u) = ^^^^ ^ + 0{l/u), 

one can express Lz{u) in terms of m^. The z-components of classical spins, are residues 
of Lziu) at u = ej (see Eq. ()3.2|) ) 



^3 = 



fefi- ^ {ej - Uk)L'_{uk) 



(4.11) 
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It also follows from Eq. ()4.2|) that f| = L^('Uj). This allows us to express variables Vj 



through Uj and the integrals of motion Hj 



{2uj - AHn 



9" 



+ 



^9 



n-l 

+ E 

fc=0 



2Hk 



g'^iuj - ek) (uj - ek)^ 



j = 0,...,n-l 



(4.12) 



(4.13) 



Thus, to determine the evolution of Sj{t) = {Kj{t)) and b{t) = {b{t)) we only need to derive 
and solve the equations of motion for Uj{t). 



V. EQUATIONS OF MOTION FOR SEPARATION VARIABLES 

In order to derive the equations of motion for separation variables Uj, we first determine 
the brackets ui^k = {Hk, ui} for /, = 0, . . . , ra — 1 using Eqs. ()4.12|) and ()4.4|) . This is done by 
taking Poisson brackets of both sides of Eq. ()4.12|) with ui and solving the resulting system 
of linear equations for ui^k- 

As soon as ui^k are found in terms of Uj and Vj, we can use the expression ()3.7j) for the 
Hamiltonian in terms of Hj to determine uj 



Ui 



{H,Uj} = J2uj,k- (5.1) 



k 



Taking Poisson brackets of both sides of Eqs. ()4.12p and ()4.13|) with ui, we obtain 

- ig'^Vj5ji = J2 — jj = 0,...n-l. (5.2) 
In order to determine Un,k and ui^k for I, k = 0, . . . ,n — 1, we need to invert the matrix 

Ajk = (5.3) 

Uj — ek 

which is called the Cauchy matrix. 

This can be done with the help of the following identity coming from partial fraction 
decomposition of the LHS: 

U?=o\u-ui) g 1 Ui^^{e,-ui) ^ ^^^^^ 



{u - Up) U?=o{u -ei) h)U- ej UiM^j - q) ' 
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To verify this identity note that both sides have the same poles and that the residues at 
these poles coincide. Substituting u = Ug in Eq. ()5.4|) . one derives 

5,, = E ' niL.K-^^) n (5.5) 

Uq - ej rii^ilei - q) f^^ Up - m 

Hence 

^ UiM^j -ei) ^_,luk-ui 

It follows from Eqs. (Q and (jO)) that 

ri/^jlMi 2&(efc -Mj) 
Finally, noticing that = L_{uj) = 2b/ g + J2k ^kl i.'^j ~ ^k) and using Eqs. ()5.ip and ()4.12p . 
we derive 



= "frT — I V i = 0> • • • > ^ - 1 



^ n—1 n—1 

b = -2ib I - + E - E "fe 

^ j=0 k=0 



(5.8) 



where the spectral polynomial Q2n+2iu) is defined as 



n-l 



Q2n+2{u) = /L2(n) n - e,)'. (5.9) 

j=0 

By Eq. (|3.4|) . the coefficients of Q2n+2iu) depend only on the integrals of motion Hj and 
parameters and g. The equation of motion for m„ = 6 is obtained by substituting Eq. ()4.10|) 
into Eq. (lO) . 

An almost identical derivation of equations of motion for separation variables can be 
performed for the mean-field BCS model (j3.1(J|) using its vector Lax representation (j3.11|) 
leading to the following equations of motion: 



2Vgff/2(^.) . , 

Uj = - ' r J=0, 

1 Im^j ) 

(5.10) 

(n n—1 ' 

gJz + E ~ E 
jr=0 k=0 



where the spectral polynomial for the BCS model is defined as 



1 n 

9 j=o 
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VI. THE SOLUTION 



Here we obtain the general solution for the mean-field dynamics of the fermion-oscillator 
(spin-oscillator) model (jl.lll.3j) by mapping it onto the corresponding BCS model. The 
solution therefore can be read off from the known general solution of the mean-field BCS 
problem ^j. 

By comparing equations of motion ()5.8|) for the spin-oscillator model with those for BCS 
model ()5.10|) . we observe that they coincide upon a replacement 

n n-1 ^ 

Y.^j+9Jz^Y.^^+n^ -J-^b, Q^c;s^(u)-^Q2n+2iu). (6.1) 

Thus, the mean-field dynamics of the spin-oscillator model p.3j) with n spins and a 
bosonic field coincides with that of the BCS model ()3.9|) with n + 1 spins. This allows us to 
immediately write down the solution for the time-dependent averages Sj(t) = (Kjit)) and 

b{t) = m) 

Ki) = [J-]bcs^ [Sj(i)]Dicke = [Sj(^)]ijc5- (6-2) 

The explicit form of [J-]bcs [^ji't)]BPS in terms of hyperelliptic functions can be found 
in Ref. ^ (see Eqs. (3.22 - 3.24) of Ref. fl). 

The dynamics of Sj{t) = {Kj(t)) and b{t) = {b(t)) is typical of a classical integrable system 
with n+1 degrees of freedomj36| (recall that n is the number of non-degenerate single particle 
levels ej in the fermion-oscillator model (jl.ll) ). The typical motion is quasi-periodic with 
n + 1 incommensurate periods. For example, a Fourier transform of b{t) = (b{t)) shows n + 1 
basic frequencies. The system uniformly (ergodically) explores the invariant torus - the 
n + 1-dimensional portion of the 2n + 2-dimensional phase space allowed by the integrals of 
motion ()3.5|) and ()3.6p . returning arbitrarily close to the initial point at irregular intervals. 

We found that in the thermodynamic limit n —>■ oo the solution simplifies for most 
physical initial conditions. The return time diverges in this limit. The dynamics of b{t) 
is particularly simple - it decays, typically as a power-law, to a steady state where \b{t)\ 
is either constant or is characterized by a few independent frequencies 311]. However, the 
motion of the spin system still contains a continuum of frequencies. The final steady state 
of \b{t)\ depends on the initial conditions. 
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VII. FEW SPIN SOLUTIONS 



The evolution described in the previous section occurs for most initial conditions and is 
stable against small perturbations of the Hamiltonian by the KAM theorem even if these per- 
turbations destroy the integrability. However, there always exists a set of points of measure 
zero in the phase space, where the motion is characterized by only a few incommensurate 
frequencies, while the stability is not guaranteed 36] . Below we consider a family of such 
solutions, which we call m-spin solutions with m < n. The reason is that in these cases the 
dynamics of n spins and the bosonic field in Hamiltonian ()1.3|) degenerates to that of m < n 
spins and the bosonic field. 

Few spin solutions are constructed by choosing integrals of motion Hj (i.e. the initial 
conditions) so that 2{n — m) roots of the spectral polynomial Q?to+?(^) defined in Eq. ()3.4|l 
become double degenerate. Suppose u = Eq is a double rootj43[ of Q2n+2iu) and note that, 
since Q2n+2(-Eo) = L^(-Eo) = 0, setting e.g. u„_i(t) = Eq solves the equation of motion 
(I5.8|l for the variable Un-i- Therefore, we can "freeze" one of the separation variables in this 
root. Then L_{Eq) = 0, and it follows from L^(u) + L_{u)L+{u) = L^{u) that L^{Eo) = 0. 
Thus, one can factor out (m — Eq) from all components of the Lax vector and show that it 
is proportional to the Lax vector of the spin-oscillator model with n — 1 spins. Similarly, if 
there are 2 pairs of double degenerate roots, the number of spins reduces to n — 2 etc. This 
procedure is followed in detail in Ref. using a different method. 

First, let us consider the general case when the number of spins effectively reduces from 
n to m < n. We have 

L(u)= (l + E^|Lt(u), (7.1) 

where dj are time- independent constants and Lt('u) is the Lax vector of the m-spin problem 

/ 2gb, \ 

2gby 
\^2u — uj' J 

The m-spin system has its own m arbitrary "energy levels" /i^. Its dynamics is governed by 
the spin-oscillator Hamiltonian ()2.3p for m spins with new parameters replacing e^- and uj. 

m— 1 m—1 

Hm=Y. + io'bb + (7 E {btk + btt) ■ (7.3) 

k=0 k=0 
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Lt(n) = 4 
9^ 



m—1 J. 

+ E (7-2) 

t^QU-^ik 



Matching the residues at poles at m = ej on both sides of Eq. ()7.ip . we express the original 
spins Sj in terms of collective spins 



(7.4) 



Constants dj are determined from the above equation using = 1/4. We have 



±1. 



(7.5) 



Finally we have to match the residues aX u = jik and the u ^ oo asymptotic. This leads 
to the following m + 1 equations 



1 + E 



j=0 l^k 



fc = 0, ...,m-l 



(7.6) 



n-1 



U = 00 — 



3=0 



Equations ()7.6|) constrain the lengths of new spins t^, i.e. the coefficients of the spectral 
polynomial Q2m+2{u) = g^hKu) Y\k=Q (u — fik^ of the m-spin system (see Eq. ()3.4p ). Indeed, 
using Eq. (j7.5|) . one can cast the constraints Eq (j7.6|) into the form 



To obtain m-spin solutions explicitly, one has to choose parameters /i^, resolve the m + 1 
constraints ()7.7j) and ()7.8|1 for the lengths of spins and frequency u', and solve for the 
dynamics of the m-spin Hamiltonian ()7.H|1 . The dynamics can be obtained from the general 
solution ()6.2j) by replacing n —>■ m and the set of {ej} with {yUfc}. 

Let us illustrate the construction of m-spin solutions by considering the cases m = and 
m = 1 in more detail. 

m=0. The 0-spin solutions correspond to the mean-field eigenstates of quantum Hamil- 
tonians (jl.lll and They are equilibrium states for the classical Hamiltonian ()2.3|1 . 

For m = the dynamics is governed by the Hamiltonian Ho = uj'hh with an obvious 
solution h = feoc"*'-'^'*''''^-', where 6o = |^(^)| = const. Further, using Eqs. ()7.2|) . ()7.4|) . and 




—6ira / = l,...,m 
9^ 



(7.7) 




(7.8) 
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()7.5|) . we obtain 

LtH = ^(4« + (2«-^')'), (7.9) 

c,iC,T = s- = , {nj) - 1 = 2s' = ^ > . 7.10 

' ^Whl + (26, - u'Y ' ^Ag%l + {2e,-uj'y 

The frequency uj' and the gap Aq = gho have to satisfy a single constraint (j7.8p that now 

reads 

^'=^-i: i (7.11) 



,.0 Jig%l + (2€, 



'^2 



We see from Eq. fl7.10|) . that all spins rotate uniformly around the 2;-axis with a frequency 
uj' . This rotation can be eliminated, i.e. uj' can be set to zero, by an appropriate choice 
of the chemical potential. Then, Eq. (j7.1ip is the analog of the BCS gap equation. The 
configuration of spins (j7.1(J|) corresponds to the product BCS wavefunction. The latter can 
be straightforwardly reconstructed from the knowledge of Sj = (K,) as fixing the average of 
spin-1/2 uniquely fixes its quantum state. The choice of signs = +1 for all j corresponds 
to the ground state, while choosing one of the signs to be —1 corresponds to an excited state 
- pair excitation of the fermionic condensate. 

There is also an important class of 0-spin solutions that is obtained by setting 6o = 0. We 
see from Eq. (j7.1(jp that in this case all spins are along the 2;-axis, s| = ej/2 = ±1/2. The 
constraint ()7.1H) is now irrelevant, because all xy-components vanish. The choice of signs 
Cj = sign ej, where are counted from the Fermi level, yields the Fermi ground state. Other 
choices correspond to the excitations of the Fermi gas. All these states are stationary with 
respect to the mean-field dynamics. For a finite number of pairs, they are non- stationary 
with respect to the quantum Hamiltonian (fTT ]) . Thus, their initial evolution is entirely 
governed by quantum corrections (cf. Refs. Iioi 111. Is^. 

m=l. In this case the dynamics reduces to that of a single collective spin t coupled to 
an oscillator, i.e. it is governed by the classical counterpart of the Dicke model ()1.4|) 

Hi = 2nt^ + uj'bb + g{bt^ + bt+). (7.12) 

Using Eq. ()7.10|) . one can express original spins in terms of the collective spin and the bosonic 
field 

^_ ^ ejjej - fi)b + gejt^ ^ ej{2ej - uj'){ej - fi) + gH, 
2gjQ^) ' 2g^jQ^ 
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To complete the construction of 1-spin solutions, we need to choose a positively defined 
polynomial Q^iu) so that it satisfies two constraints ()7.7p ()7.8|) and solve for the dynamics 
of the classical Dicke model (|7.12p . In this case there is only one non-stationary separation 
variable Uq and equations of motion (j5.8j) take the following form: 



+ = 0, 

b = —ih{uj + 2/i — 2uq). 



(7.14) 



These equations can be solved in terms of elliptic functions. This is not surprising, since 
the model ()7.12p is, in fact, equivalent to a spherical pendulum and its solution has already 



been obtained in Refs. 



3, 



34|. 



where they were used to 



The 1-spin solutions were originally obtained in Refs. 
describe the evolution beginning from a state infinitesimally close to the normal ground 
state. 



VIII. CONCLUSION 

In this paper we solved for the mean-field dynamics of the fermion-oscillator model (jl.lj] . 
In the mean-field approximation, the problem reduces to a classical Hamiltonian model. We 
derived integrals of motion for both the classical ()2.3|) and the quantum (ll.lj) models and 
showed that the dynamics of the fermion-oscillator model (jl.ll) maps onto that of the BCS 
model. This was used to derive an explicit general solution for the mean-field dynamics of 
the fermion-oscillator model with an arbitrary finite number n of degrees of freedom. 

The typical dynamics is quasi-periodic with n incommensurate frequencies. The system 
ergodically explores the part of the phase-space allowed by integrals of motion, returning 
arbitrarily close to its initial state at irregular time intervals. We have also constructed a 
class of particular, few-spin, solutions, for which the dynamics reduces to that of m collective 
spins governed by the same classical model. The case m = corresponds to the mean-field 
eigenstates of the fermion-oscillator Hamiltonian (jl.lj) . 

An interesting problem is to obtain and describe the thermodynamic limit, n — 00, of 
the solution. For most physical initial conditions the dynamics considerably simplifies in 
this limit [3^. In particular, we found that decays, typically as a power-law, to a 

steady state where \(b{t))\ is either a constant (cf. Refs. 12, y, LD, 111! ) or oscillates with few 
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independent frequencies (cf. Fig. 2 in Refs.y,|6|). The motion of the spin system still contains 
a continuum of frequencies. The steady state of \{b(t))\ depends on the initial conditions. 

It is also interesting to study quantum and finite temperature effects. Quantum correc- 
tions to the mean-field dynamics can become important when the bosonic mode is weakly 
populated 3^. For example, normal states where (b'^b) = {cjiCj^) = are stationary in the 
mean- field approximation ()2.2|) . but are not stationary with respect to quantum dynamics 
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